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A BCS CONDENSATE IN NJL 3 + i? 
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We present results from a lattice Monte Carlo study of the Nambu - Jona-Lasinio 
model in 3+1 dimensions with a baryon chemical potential p A 0. As p is increased 
there is a transition from a chirally-broken phase to relativistic quark matter, 
in which baryon number symmetry appears spontaneously broken by a diquark 
condensate at the Fermi surface, implying a superfluid ground state. Finite volume 
corrections to this relativistic BCS scenario, however, are anomalously large. 


The proposal that at high baryon density a condensation of diquark 
pairs in degenerate quark matter takes place, leading to color supercon¬ 
ducting ground states, has generated a profusion of theoretical and phe¬ 
nomenological interest 1 . Most approaches assume the existence of a Fermi 
surface and then consider its instability with respect to diquark conden¬ 
sation arising from an attractive qq interaction. The resulting excitation 
energy gap A can be calculated self-consistently in close analogy to the 
BCS gap in electronic superconductors. Model calculations yield gaps as 
large as 50 - lOOMeV 2 . 

A first-principles calculation of the ground state of nuclear or quark 
matter using the methods of lattice QCD, however, is hampered by the 
“sign problem”, ie. exp (—S/h) is no longer positive definite, and cannot 
be used as a probability measure in importance sampling the Euclidean 
path integral, once baryon chemical potential /.i ^ 0. It is still possible, 
however, to study simpler models of strong interactions with fi yf 0 using 
lattice methods 3 . The Nambu Jona-Lasino (NJL) model, which closely 
resembles the original BCS model, is particularly interesting in the current 
context, because it is the only simulable model which exhibits a Fermi 
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surface for sufficiently large /i. The Lagrangian density is 

L = ip( $ + mo + o)ip - y [W) 2 - (V>75 rip) 2 ] . (1) 

Here tp and ip are isospinor “quark” fields carrying baryon charge; note that 
both qq and qq interactions are attractive in the isoscalar channel. 

For coupling strength g 2 sufficiently large the vacuum ground state has 
chiral SU(2)L®SU(2)fl symmetry spontaneously broken to SU(2)j SOSp i n by 
generation of a chiral condensate and resulting “constituent” quark mass 
£ ~ g 2 {ipip). At zero temperature this state is stable as g is increased 
until at n c ~ £ chiral symmetry is restored, at which point the baryon 
density ng = (ip'yoip) rises from zero 4 . For /i > [i c the ground state is 
thus relativistic quark matter with Fermi momentum fcp ~ g. A natural 
question is whether qq attraction in this phase leads to a diquark condensate 
spontaneously breaking the U(l)s symmetry of (1) leading (since in this 
case the broken symmetry is global rather than local) to superfluidity. 

The NJL model in 2+1 dimensions has been studied extensively using 
Monte Carlo simulations. While there is enhanced qq pairing in the scalar 
isoscalar channel for g > qi c 5 , simulations with an explicit U(1 ^-violating 
source term jqq found no evidence for diquark condensation in the j —> 0 
limit 6 . Instead, the condensate scales approximately as ( qq(j )) oc j+, 
the exponent d increasing with g. Moreover, the quasiparticle dispersion 
relation in the spin-i channel is consistent with a vanishing BCS gap A = 
0. The interpretation is that long-wavelength fluctuations of the complex 
phase of the qq wavefunction destroy the condensate, but that long-range 
phase coherence remains (ie. correlations are dominated by “spin-waves” 
and remain algebraic functions of spatial separation) just as in the low- 
T phase of the 2d XY model. Superfluidity in two dimensions is realised 
in Kosterlitz-Thouless mode; namely, persistent flow patterns can only be 
disrupted by creating a vortex/anti-vortex excitation and then translating 
one member of the pair around the universe, which in the low-T phase 
entails an energy cost increasing logarithmically with system size. 

This rather exotic scenario, namely gapless thin film superfluidity, is in¬ 
herently non-perturbative and is not exposed by self-consistent approaches. 
It thus justifies the application of numerical simulations, and begs the 
question of whether similar surprises will emerge from a study of the 
phenomenologically-motivated case of NJL 3 + i 7 . In 3+ld the model is 
non-renormalisable, and hence its physical predictions are sensitive to the 
details of the ultra-violet cutoff. We have matched our lattice model to 



February 1, 2008 19:35 WSPC/Trim Size: 9in x 6in for Proceedings 


heidelberg 


3 

low-energy QCD following the procedure of Klevansky 4 , working to lead¬ 
ing order in an expansion in 1 /N c , where N c is the number of quark “colors”, 
and calculating E, m v and in terms of bare parameters g 2 and mo- Since 
we use staggered lattice fermions y, y, we have N c = 4, in which case rea¬ 
sonable low-energy phenomenology emerges for the choice a 2 / g 2 = 0.565, 
TOoa = 0.002, with inverse lattice spacing a -1 ~ lGeV. Figure 1 shows 



Figure 1. Chiral condensate {xx)> baryon density, and diquark condensate as functions 
of /x, showing both large- N c solution (solid curve) and simulation results (points). 


results for (xx) an d and confirms the transition between vacuum and 

2 

quark matter at g c ~ E ~ ^-(yy). The discrepancy between line and 
points is ascribed to 0(1/N c ) corrections and amounts to a 30% effect. In 
contrast to 2+1 d where a sharp first-order transition is observed 3 ’ 5 , the 
transition appears continuous in the chiral limit too —> 0, although this may 
be sensitive to the cutoff 4 . 

To investigate the diquark sector, we introduce a source term j±(q tr q± 
qq tr ) into the quark operator used for measurements (since the simulation 
dynamics are performed with j± = 0, this is a “partially quenched” ap¬ 
proximation) and define the following: 


(QQ±) = 


1 5 In Z 

v~dJ7 


d{qq±) _ 
dj± 


R(h) 


X+ 


: x± 


X- 
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The nature of the ground state is then determined by 

lim (qq+) = 0, lim R = 1 if U(1 )b manifest 

3 +^ 0 J+-* 0 

lim (qq+) 0, lim R = 0 if U(l)s broken (3) 

o j+^o 

where the vanishing of R in the broken phase follows from the Goldstone 
nature of the qq- boson in this case. 



j i 


Figure 2. Susceptibility ratio R (left) and diquark condensate (qq+) (right) as functions 
of j-\- for various [i. 


Figure 2 shows results for R and (qq+) from data taken on 12 4 , 16 4 
and 20 4 lattices and then linearly extrapolated to the limit L )" 1 = T —> 0. 
For n = 0 the data appear to extrapolate smoothly to R = 1, (qq+) = 0 
as j+ —> 0. For fia = 1.0 by contrast, which lies in the region of high 
baryon density and restored chiral symmetry, then if data with j + a < 0.2 
are ignored the extrapolated values are R ~ 0, ( qq+)a 3 ~ 0.28, suggesting 
superfluidity. For intermediate /i the situation is not so clear cut, but 
the results of a quadratic extrapolation of (qq+ ) to = 0 are shown in 
Figure 1. For 0.4 < Ra 0.9 the condensate scales as r 2 , which is expected 
if diquark pairs within a distance A of a Fermi surface of radius r participate 
in the condensate leading to (qq+) ~ A . The value of (qq+) at jia = 1.0, 
together with an estimated E ~ 350MeV, is thus plausibly consistent with 
the gap values reported in [2]. The positive curvature contrasts with the 
behaviour of the superfluid condensate in Two Color QCD in which qq pairs 
are tightly bound and form a conventional Bose-Einstein condensate 8 . 

Since our interpretation of broken U(l )b symmetry at high baryon den¬ 
sity relies on discarding some data, the above conclusions are necessarily 
provisional. As spontaneous symmetry breaking does not occur in a finite 
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volume, it is natural to attribute the depletion of (qq+) at small j + to the 
effects of a finite L s . Note however that the threshold ja ~ 0.3 required for 
behaviour characteristic of the thermodynamic limit seems anomalously 
large compared with that required to observe chiral symmetry breaking 
(moa ~ 0.002) on the same volume, or indeed for thin film superfluidity 
(ja ~ 0.1) in NJL 2 +i 6 . We suspect that finite volume effects in systems 
exhibiting a Fermi surface are unconventional, and indeed data on systems 
with varying L s ,L t suggest the approach to the thermodynamic limit at 
small j is non-monotonic '. Solution of the gap equation on finite systems 
shows that A (L s ) oscillates before approaching its infinite volume limit for 
L s <1 8fm, which translates to a 40 3 spatial lattice for our model 9 . 

In future work we plan to study the quasiparticle spectrum which in 
principle permits a direct estimate of the superfluid gap A. It might also 
prove interesting to choose different bare parameters, since a small change in 
the cutoff might actually change the chiral symmetry restoring transition to 
first order 4 , making it easier to distinguish the properties of high- and low-/i 
phases. Finally, if our conclusion that NJL 3 + i is superfluid at high density 
perists, then it will be interesting to examine the stability of the condensate 
as either temperature or fj,i SO spin (which has the effect of separating u and 
d Fermi surfaces and hence inhibiting pairing in the isoscalar channel) is 
raised from zero. 
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